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' Abstract. The continuum and emission-line spectrum of the narrow-line Seyfert 1 galaxy 

M-T 

Ark 564 is used to investigate, for the first time, the physical conditions and structure of its 

O ; 

narrow line region (NLR). For this purpose, composite models, accounting for the coupled 
effect of photoionization and shocks, are employed. The emission-line spectrum of Ark 564, 
which ranges from the ultraviolet to the near-infrared, shows a rich forbidden line spectrum. 
Strong emphasis is given to the study of the coronal line region. The diversity of physical 
conditions deduced from the observations requires multi-cloud models to reproduce the ob- 
served lines and continuum. We find that a combination of high velocity (V,= 1500 km s -1 ) 
shock-dominated clouds as well as low velocity ( V s = 150 kms~') radiation-dominated 
clouds explains the coronal lines, while the optical low-ionization lines are mainly explained 
by shock-dominated clouds. The results for Ark 564 are compared with those obtained for 
other Seyfert galaxies previously analyzed such as NGC5252, Circinus, NGC4051 and 
NGC4151. The model results for the ultraviolet and optical permitted lines suggest that 
the broad line region may contribute up to 80%, depending on the emission-line, being of 
about 30% for Hy8. The consistency of the multi-cloud model is checked by comparing the 
predicted and observed continuum, from radio to X-ray, and indicate that the dust-to-gas 
ratio in the clouds varies from 10~ 15 to 10~ 12 . 



Key words, galaxies: Seyfert - galaxies:individual: Ark 564 - infrared: galaxies - 
line:formation 



2 Marcella Contini et al.: The Narrow Line Region of Ark 564 

1. Introduction 



Narrow-line Seyfert 1 galaxies (NLS1) are a peculiar group of active galactic nuclei (AGN) with 
extreme properties. Their permitted optical lines emitted by the broad line region (BLR) show 
full widths at half-maximum (FWHM) not exceeding 2000 kms -1 , the [Oiii]fH/3 ratio is less 
than 3, and the UV-optical spectrum is usually dominated by strong Fen emission multiplets 
(Osterbrock & Pogge lT985> . In the soft X-ray, NLS1 have generally much steeper continuum 
slopes and rapid variability (Boiler et al. 1996 1. In addition, Leighly ( 1999 1 found that the hard 
X-ray photon index is significantly steeper in NLS 1 compared with that of classical Seyfert 1 
and that soft excess emission appears considerably more frequently in the former than in the 
later. Nowadays, the most accepted explanation for the above properties is that these objects 
have relatively low black-hole masses for their luminosity but high accretion rates (Boiler et al. 
I1996> . In fact, observational evidence obtained through reverberation mapping techniques points 
toward the presence of a black hole of smaller mass in NLS1 (Mbh ~ 10 6 ~ 7 M ) than that of 
other broad-line AGN (M B h ~ 10 7 ~ 9 M ; Peterson & Wandel 2000 1 with similar luminosity. 

While most efforts have been devoted to explain the extreme properties associated to the BLR 
in NLS1, very little has been said about the properties of the narrow line region (NLR) in these 
objects and how do they compare to that of other AGN, except for a few works reported in the 
literature (Goodrich [19891 Rodriguez- Ardila et al. 120001 Nagao et al. 120011 Rodriguez- Ardila 
et al. I2002> . Given that the NLR is one of the fundamental ingredients of AGN, every effort 
to understand what NLS 1 are must include a detailed investigation of the NLR at least for the 
following three reasons. 

First, the NLR is the largest structure directly affected by the radiation from the central en- 
gine. Therefore, if the spectral energy distribution (SED) of the ionizing radiation flux is some- 
how different from that of classical Seyferts, the resulting ionization structure of the NLR gas 
should be different, as it is observed in some NLS 1 (Rodriguez- Ardila et al. l2000l 

Second, NLS 1 are located at the extreme end of the anticorrelation between [O m] and Fe n 
(Boroson & Green 1992 ), namely, Fe n is strong in objects with weak [O m] emission and weak in 
those with strong [O m]. That inverse correlation suggests the existence of a physical mechanism 
which balances Fe n excitation against the illumination of the NLR. Although the nature of the 
Fe n is out of the scope of this paper, the study of the NLR in NLS 1 may provide important clues 
to identify the parameter linking BLR and NLR excitation. 

Third, the spectra of NLS 1 usually display strong high ionization forbidden emission lines 
(also known as coronal lines) such as [Fex] /16374, [Fexi] /17892 and [Fexiv] /15303 (Penston 
et al.[l984; De Robertis & Osterbrock 1984). Up to now, it is not clear what the physical mech- 
anisms involved in the production of these features are, but three possible scenarios have been 
suggested: (i) collisional ionization of gas with temperatures of T e ~ 10 6 K (Oliva et al. 119941 : 
( ii) photoionization by the central non-thermal continuum (Ferguson et al. U997t : and ( Hi) a com- 
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bination of shocks and photoionization (Viegas-Aldrovandi & Contini 1989, Contini & Viegas 
2001 1. Although these lines are common in most AGN, NLS1 appears as suitable candidates to 
investigate their origin because they are not severely blended with nearby broad features as in 
classical Seyfert 1 or diluted by stellar absorption features as in Seyfert 2. 

The study of the NLR in NLS 1 is not straightforward, however. Deblending the optical per- 
mitted lines in NLS 1 is difficult because no transition between the narrow and broad components 
is observed. This shortcoming translates into large uncertainties in determining the fraction of 
Ha and H/3 that is emitted by the NLR, affecting in turn, the determination of the reddening of 
the NLR gas. The extinction correction should then relay on other indicators such as ratios of for- 
bidden lines of the same ion, separated by a large wavelength interval, yet not always observed 
with the same instrumental setup. 

For the above reasons, we have started a program to study, in a consistent way, the NLR 
of a selected sample of NLS1 galaxies with the aim of investigating the physical conditions 
and structure of that region, with particular emphasis in the origin of the coronal lines. To that 
purpose, composite models accounting for both the effects of photoionization by a central con- 
tinuum source and shocks, created by collision of gas clouds with the intercloud medium, will 
be employed. The well-known NLS1 galaxy Arakelian564 (z = 0.0247) is a suitable candidate 
to begin with. The strong coronal line spectrum observed by Rodriguez- Ardila et al. (2002] here- 
after RVPP02) in the near-infrared region (NIR) as well as the optical measurements reported 
by Erkens et al. ( 1997 hereafter E97) on this source form one of the most complete set of lines 
available for modeling the NLR and studying the effects of shocks on line intensities in a NLS 1 . 
Following Moran (2000), Ark 564 exhibits an unresolved radio core and what appears to be a 
small scale (~ 1") jet. This feature strengthens the hypothesis that shocks are at work. Moreover, 
it suggests that extended synchrotron radiation should be observed in the radio range. 

Ark 564 has been extensively studied from radio to X-rays, allowing a consistent modeling 
of both the line and continuum emission. Among other reasons that turns this target interesting 
is that it is the brightest known NLS1 in the 2-10 keV band (Collier et al. l2001i . The results of 
an intensive variability campaign on several wavelength bands show that the optical continuum 
is not significantly correlated with the X-rays (Shemmer et al. 2001 1. The ionization state of the 
gas, as described by Crenshaw et al. d2002t . is relatively high. They found that at least 85% of 
the narrow emission-line flux comes from a region < 95 pc from the nucleus and surrounded by 
a dust screen associated to a "lukewarm" absorber. 

In Sect. 2 we describe the data set used to the study the NLR. A description of the modeling 
output and fitting to the line fluxes is given in Sect. 3. As a consistency check, the fit of the 
continuum emission SED is presented in Sect. 4. Discussion of results follows in Sect. 5 and 
Conclusions in Sect. 6. 
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2. The data 

The NLR emission-line fluxes used in this work were primarily taken from E97 and RVPP02 
for the optical (0.34-0.95 fim) and NIR (0.80-2. 34//m) lines, respectively. Both sets of data were 
taken at similar resolution and S/N. Although the slit width employed in the optical observations 
was larger than that of the NIR ones (1.5" and 0.8", respectively), the contamination introduced 
by the host galaxy or extended emission to the optical fluxes is negligible. In both cases, the 
authors report that most of the integrated signal came from the unresolved nuclear region and 
that the spectra were extracted from pixels centered on the maximum of the continuum emission. 
This is confirmed by comparing the flux of the [S ra] 0.953 yum line, common in both spectra. 
E97 measured a value of 3.06+0.13 x 10~ 14 erg s _1 cirT 2 , in very good agreement with the value 
of 3. 1 1 ±0. 10 x 10" 14 erg s" 1 cnr 2 reported by RVPP02. 

For reddening correction, we used the value of E(B-V)=0. 09+0. 03 reported by E97. Although 
it is a bit lower than the E(B-V)=0. 14+0.04 reported by Crenshaw et al. (2002), they agree 
within errors. The former value was determined applying the method of Allen ( 1979) to the [S n] 
0.4068yum, 0.4076/mi, [Sn] 0.6717^, 0.6731//m, [On] 0.3727 /mi and [On] 0.7330/mi line 
pairs, and so, most suitable to the NLR. Cols. 2 and 5 of Table[T]list the extinction corrected line 
ratios relative to the total H/3 flux. In order to avoid the uncertainties introduced by the deblending 
of the observed permitted lines into the component emitted by the NLR and BLR, we decided 
to normalize the narrow line fluxes to that of [Om] 0.5007+0. 4959/mi. The resulting values are 
reported in Cols. 3 and 6 of that same Table. 

3. Modeling the spectra 

The main goal of this work is to model the NLR spectrum of the NLS1 galaxy Ark 564 using 
the fluxes of optical and NIR coronal lines (see Tabled as major constraints. This ensures that 
important information about high ionization material (see Contini et al. 2002a) is taken into 
account. The data from Table ^ imply that gas clouds with very different physical conditions 
must coexist in the NLR of this object. This is reflected, e.g. in the large interval of critical 
densities (10 3 — 10 9 cirT 3 ) of the transitions and ionization potentials (up to 370eV) of the ions 
involved in the production of the observed forbidden lines, or in the large range of line widths 
-a factor of 3 in FWHM (see Table 5 of RVPP02) measured in the line profiles. With this in 
mind, we have run a grid of models with parameters in the range suitable to the NLR of AGN 
and selected the models that best matches at least some of the observed line ratios. The grid 
was based on Contini & Viegas (2001 ). The results are obtained by the SUMA code (Viegas & 
Contini H"994> which, besides the line and continuum spectra emitted by the gas, also calculates 
re-radiation by dust. 

Several multi-cloud models, obtained from a weighted average of single clouds models, were 
calculated and compared with the observed data. For each multi-cloud model, higher weights 
were assigned to single-cloud models matching the largest number of line ratios. Using this 
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Table 1. Emission line ratios 1 for Ark 564, corrected for E(B-V)=0.09 
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Line 




^/[Oin] 


Line 




F,/[Om] 


[Ne v] ,13425 


16.5 ±0.9 


13.9 ±0.9 


[Fex] ,16374 


7.10 ±0.42 


5.97 ± 0.41 


[O II] ,13727 


16.3 ±0.7 


13.7 ±0.7 


[Nil] ,16548 


0.673 ± 0.041 


2.79 ±0.19 


[Fev] ,13839 


0.226 ± 0.028 


0.19 ±0.02 


Ha 11 


355.0 ± 14.7 


298 ± 16 


[Nem] ,13868 


9.75 ± 0.41 


8.2 ±0.4 


[Nil] ,16583 


9.97 ± 0.44 


8.30 ± 0.46 


[S ii] ,14068 


1.96 ±0.11 


1.65 ±0.11 


He i ,16678 s 


1.73 ±0.13 


1.46 ±0.12 


[S ii] ,14076 


0.482 ± 0.038 


0.40 ± 0.03 


[Sn] ,16716 


3.88 ±0.32 


3.26 ± 0.29 


[0 in] ,14363 


4.89 ±0.21 


4.11 ±0.22 


[S ii] ,16730 


4.51 ±0.20 


3.79 ±0.21 


He II ,14686" 


14.8 ± 4.7 


12.4 ± 3.9 


[Arv] ,17005 


0.656 ± 0.061 


0.55 ± 0.05 


[Anv] ,14711 


1.01 ±0.09 


0.85 ± 0.08 


[Arm] ,17135 


1.69 ± 0.09 


1.42 ±0.09 




100.0 


84.1 


[Arm] ,17751 


0.538 ±0.100 


0.45 ± 0.09 


[Fevn] ,14893 


0.573 ± 0.075 


0.48 ± 0.06 


[Fexi] ,17891 


8.31 ±0.39 


7.0 ± 0.40 


[Fevil] ,14988 


0.573 ± 0.088 


0.48 ± 0.08 


[S m] ,19069 


5.76 ± 0.26 


4.80 ± 0.27 


[Om] ,15007 + 4959 


118.9 ±4.0 


100.0 


[S m] ,19532 


10.7 ±0.5 


9.00 ± 0.52 


[Fevn] ,15145 


0.283 ± 0.023 


0.23 ± 0.02 


[C I] 0.985 pm 


0.23 ± 0.07 


0.19 ±0.06 


[Fevi] ,15159 


1.01 ±0.09 


0.85 ± 0.08 


[S vm] 0.991 pm 


1.75 ±0.13 


1.47 ±0.12 


[Ni] ,15199 


2.04 ±0.12 


1.72 ±0.12 


[Fexm] 1.074 pm 


2.88 ± 0.55 


2.42 ± 0.47 


[Fevi] ,15335 


0.697 ± 0.097 


0.59 ± 0.08 


[Six] 1.252 pm 


1.99 ±0.19 


1.60 ±0.16 


[Fevn] ,15721 


2.31 ±0.15 


1.94 ±0.14 


[Fen] 1.257 pm 


0.97 ±0.16 


0.82 ±0.14 


He i,15875 a 


9.92 ± 0.52 


8.34 ±0.52 


[Six] 1.430 pm 


4.55 ± 0.32 


3.83 ±0.30 


[Fevn] ,16087 


3.64 ±0.19 


3.06 ±0.19 


[Si VI] 1.963 pm 


1.89 ±0.35 


1.59 ±0.30 


[Oi] ,16300 


5.29 ± 0.23 


4.45 ± 0.24 


[Cavm] 2.321 pm 


0.82 ± 0.26 


0.69 ± 0.22 



1 Fluxes in the optical region, up to 8000 A, were taken from E97. NIR lines, from 9000A and longward 

were taken from RVPP02. 
" Permitted line with contribution from the BLR. 



approach, we selected two final models that emerged as the best solutions. In the following 
sections they will be analyzed in detail to determine which one best represents the NLR state of 
Ark 564. 

3.1. The models 

The general model assumes gaseous, dusty clouds, which move outwards from the active center 
(AC) and emit the line and continuum spectra. A shock front forms on the outer edge of the 
clouds, while the inner edge is reached by the radiation from the AC. 

The input parameters are: the radiation flux intensity, F H (in photons crrT 2 s _1 eV~' at 1 
Ryd), the spectral index, a (= - 1 .5), the shock velocity, V„ the preshock density, no, the preshock 
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magnetic field, Bo (= 10~ 4 gauss), the dust-to-gas ratio by number, dig, and the geometrical 
thickness, D, of the emitting clouds. Cosmic relative abundances (Allen 1973) are adopted. In 
radiation dominated (RD) clouds the effect of the flux from the AC prevails on the shock. In 
shock dominated (SD) clouds the flux is absent (F h= 0). 

Notice that photoionization models with spectral index a=-2 were adopted by Rodriguez - 
Ardila et al. ( 2000 1 in order to explain the low [O m]/HB line ratio observed in the NLS 1 galax- 
ies. The observations required both matter-bounded and ionization-bound clouds. However no 
attempt was made to fit coronal lines or other emission lines in wavelength ranges other than 
the optical one. The presence of a non-negligible velocity field, indicated by the FWHM of the 
emission-lines, suggests that composite models, accounting for shock and photoionization, and 
also providing regions of high and low density gas, are probably more appropriate to model the 
emission-line spectra of NLS1 galaxies. 

Investigating the optical and UV properties of NLS 1 galaxies, Constantin & Shields ( 2003 ) 
suggest a spectral index of -1.34 in the UV-blue continuum. On the other hand, ROSAT data 
show that the X-rays spectrum shows a power index steep and equal to -2.4 (Boiler, Brandt & 
Fink 1996 1, while ASCA data, corresponding to a different energy range and used to analyze the 
X-rays variability, indicate a power index of the order of -1.3. Analyzing the optical spectrum of 
Ark 564, Crenshaw et al. (2002) assume that the continuum of Ark 564 consists of power-laws 
with index -1.0 for energies less than 13.6 eV, -1.07 for 13.6 eV < E < 0.5 keV, and -1.6 for E > 
0.5 keV (consistent with Turner et al. 2001 ). Therefore, at least from the observational point of 
view, there is no strong reason to adopt a power index steeper than -1.6. 

Here we adopt a UV - X-rays power index a = -1.5 for all the models, recalling that the 
shocked zone also contributes to the emission-line intensities, so that our results are less depen- 
dent on the shape of the ionizing radiation. As shown in the next sections, our models with a = 
-1.5 provide a good fit to the observed emission-line spectrum in a large range of wavelengths. 

In order to account for the constraints imposed by the observations (line widths, densities, 
etc.), the models selected from the grid show a large range of shock velocities (V s = 100 - 1500 
kms -1 ) and preshock densities («o = 100 - 700 cm 4 ), as in shown in Table[2] Notice that the 
shock velocities in the models lead to temperatures in the post shock region between 1.5 xlO 5 
K (Vj=100 kms -1 ) and 3xl0 7 K (V s =1500 kms" 1 ), in agreement with the average temperature 
derived by E97 from the [Fe vn] lines for the region producing the coronal lines. An overview 
of the values for the initial conditions adopted for each of the 14 individual cloud models are 
listed in the first six entries of Table|3] Notice that the values wr correspond to the model relative 
weight. Thus, in order to obtain the fraction of the HB luminosity emitted by a given cloud, 
the relative weight should be multiplied by the corresponding HB value listed in Table |2] and 
normalized by the sum over all clouds. 

The cloud geometrical thickness D was chosen between 3 and 10 pc, but some models repre- 
sent filaments with D < 1 pc. These filaments are created by fragmentation of matter in a turbulent 
regime in the presence of shocks. Such clouds must be present in order to emit strong [Fe xi] con- 
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sistent with weak [Ci] lines. The clouds emitting strong Fe coronal lines are matter-bounded, so 
they are mostly composed of highly ionized gas at high temperature. In a search for variations of 
[Fevn] and [Fex] in high ionization Seyfert galaxies, Veilleux ( 1988) claims that neither [Fex] 
nor [Fe vn] in Ark564 show any significant variations between 1977 and 1988. Considering that 
time variations are related to D/V s (the time in which the shock crosses a filament), a minimum 
time of 25 years would be needed to detect any variability on the line flux emitted by clouds of 
model ml4 (see Table 2). This time would be much higher for other clouds (see the values of 
D/V s ), in perfect agreement with Veilleux's results. 

The flux intensities range from in SD models to Fh— 5x10 12 photons cirT 2 s eV~' at 1 
Rydin RD models. The luminosity of Ark 564 at 1 Ryd, L=7xl0 43 erg s" 1 (Romano et al. [2002l . 
constraints the models. Actually, the luminosity at 1 Ryd corresponding to the maximum F H = 
5xl0 12 photons cirT 2 s _1 eV~' is 4.5xl0 44 erg s , assuming full coverage of the source. 

The main criterion adopted for selecting the final multi-cloud models is that the ratio of cal- 
culated results to the data is within a factor of ~ 2. However, given that the clouds are distributed 
over a large volume of space, subjected to a variety of physical conditions, some line ratios may 
show slightly larger discrepancies. The single-cloud models which lead to the best two fitting 
multi-cloud models are presented in Tableland labeled Ml and M2. The resulting line ratios 
normalized to [O m] /1/15007+4959 are compared with the observed ratios in Table[3] 

The multi-cloud model Ml shows that [Ora] is overpredicted by a factor of ~ 1.5 (compare 
the resulting ratio for Ha) and [Fex]/[Fexi] is higher than observed by a factor of >2. Notice, 
however, that the observed lines (E97, RVPP02) show very complex line profiles, which could be 
an indication that more single-cloud components should be included provided they do not change 
significantly the overall spectrum emitted by the cloud combination. However, in order to include 
new components, additional constraints coming from the line profiles should be accounted for, 
which is out of the scope of this paper. Moreover, model Ml gives [Nn] and [S n] values higher 
than observed. Since N and S are not the main gas coolants, this problem could be solved by 
changing the adopted abundances of these two elements without affecting the results for other 
lines. 

On the other hand, model M2 shows a better agreement with the observed [Fe x]/[Fe xi] ratio 
(<3), even if the [Fe xi]/[0 m] ratio is underpredicted, as well as [Oi]/[Oiii], while the calculated 
He ii/[0 m] ratio is high. In addition, [N n] and [S n] ratios are too low in this model, which would 
be improved using a higher abundance for N and S. 

Both models Ml and M2 underpredict the [N i]/[0 m] ratio. Notice that the [N i] line is usually 
weak because the critical density for collisional deexcitation is very low (~ 1000 cm 4 ) and the 
volume transition zone of the clouds is usually very low or inexistent. 

Apart for the above discrepancies, it is clear from the data of Table[3]that model Ml leads to 
a better agreement than M2 for a larger number of lines. But before we can definitively choose 
either of the two models as the one that better represents the NLR state of Ark 564, further 
constrains must be applied. 
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3.2. Relative contribution to the lines from different clouds - weights and line widths 



The relative weights adopted to sum up the single-cloud models, wr, are given in Table [2] The 
relative contribution (in per cent) of each cloud to the total predicted line flux is also listed in 
Table|2]for the multi-cloud models Ml and M2. 

The physical conditions of the gas in each cloud are represented by the ensemble of the 
input parameters of the corresponding model. In composite models, the shock velocity and the 
flux from the AC are fundamental. Although the shock velocity does not correspond exactly to 
the velocity of the gas, it can be used to discuss the FWHM of the lines presented by RVPP02 in 
their Table 5. Clearly, all the line profiles show some substructure, confirming that many different 
conditions contribute to the different lines. Particularly, the FWHM are all > 300 kms" 1 . This 
is confirmed by our modeling, which shows (see Table |2j that all the coronal lines have strong 
contributions from clouds with V, > 300 km s -1 . 

As an example, the [S m] line is mainly produced by clouds with V s =500 km s -1 (50% in 
model Ml and 95% for model M2), in good agreement with the value of FWHM = 440 km s -1 
of RVPP02. The [S vm] and [S ix] lines, with more than 50% emitted by clouds with V s < 300 
kms -1 , have large wings with V,=1500km s . Regarding the [Fe vn] line, it is better reproduced 
by model M2 than by model Ml. In fact, model M2 shows that 60% of that line is emitted by 
clouds with V s =300-500 kms -1 , very close to the observed FWHM of 420 kms -1 . [Fex] and 
[Fexi] are another set of lines that deserves comments. Observationally, their line profiles have 
similar velocity distributions, with FWHM of 570 kms -1 and 695 kms -1 , respectively. This 
result must be compatible with the model outputs. In fact, model Ml predicts a large contribution 
(~60%) from high velocity clouds (V s > 500 kms -1 ) for both lines. In contrast, according to 
model M2, [Fe x] and [Fe xi] are mostly emitted by low velocity clouds with V s =100 km s -1 and 
150 km s -1 . The contribution of high velocity clouds, in this case, is not larger than 30%. 

In addition to the line broadening produced by the shock, bulk motions of the clouds around 
the gravitational potential of the central mass concentration should introduce an additional broad- 
ening, not taken into account in our modeling. Although a detailed calculation of the resulting 
profile is out of the scope of this paper, it is important to note that our results, overall, are consis- 
tent with the observed line widths. The reader should be aware that a complete modeling, able to 
simultaneously reproduce line fluxes and widths as well as local variations in most physical pa- 
rameters, is not realistic because of the large number of free parameters, some of them unknown, 
that such an approach would require. 

We conclude from the above discussion that model Ml is more representative of the actual 
conditions in the NLR of Ark 564 than model M2. In addition to better fit the observed line ratios, 
the predicted velocity distribution for the coronal gas is more consistent with the observations. 
The choice of model Ml is further reinforced by the discussion in the next section. 
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3.3. Relative contribution of the BLR emission 
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As it was pointed out in the Introduction Section, one of the main problems associated to the 
study of the NLR in NLS 1 galaxies is the difficulty in obtaining a reliable estimate of the per- 
mitted lines fluxes emitted by that region. We recall here that the models calculated by SUMA 
refer to the NLR of Ark 564 because the code is not adapted to calculations for emission from 
a very dense gas. However, comparing the observed and calculated line ratios relative to H/3, as 
is shown in Cols. 7 and 8 of Table|3]for models Ml and M2, respectively, allows us to derived 
roughly the contribution of the broad H/3. This is possible because our permitted line flux mea- 
surements include the contribution from both the NLR and BLR while the predicted line ratios, 
relative to H/3, do not include the broad line contribution. Thus, we expect that the calculated 
forbidden lines ratios, relative to H/3, be higher than the observed ones. 

This is in fact the situation found for the values predicted by model Ml. They are ~1.3 
higher relative to the ratios reported in Cols. 2 and 5 of Tabled indicating that the predicted H/3b r 
corresponds to about 70% of the total observed H/3. On the other hand, for model M2 the values 
listed in Col. 8 of Table[3]indicate the opposite, suggesting that this model is unrealistic. 

This result supports Rodriguez- Ardila et al. ( 2000 1 conclusions about the fraction of the NLR 
flux present in the observed H/3 line of NLS1 galaxies. They found that, on average, 50% of the 
total integrated H/3 flux is emitted by the NLR, in contrast to the fraction of 10% usually found 
in broad-line Seyfert 1 galaxies. 

Notice that the ratio (h/H/3) calc /(l A /H/3) obs derived for He n /14686/HjS, He i A5876/H/3, is not 
constant, indicating that the broad contribution to He n /14686, He i, and H/3 is different. 

Finally, in Table 0] we compare the results of model calculations with the UV line ratios 
observed by Crenshaw et al. (2002 1, which include both the NLR and BLR contributions. Ratios 
less (higher) than unity indicate that the contribution of the BLR to the line is high (low). 

4. The SED of the continuum 

Following the results from the previous sections, model Ml is favored and its consistency with 
the observed continuum in a large range of wavelengths is checked below. In Fig.[flthe spectral 
energy distribution of the bremsstrahlung and of re-radiation by dust from model Ml is compared 
with the observed spectral energy distribution of Ark 564. The data were taken from the De 
Vaucouleur et al. J199U : Zwicky & Kowal (19681; De Vaucouleurs & Longo ( 1988 1; Moshir et 
al. {T990l and Dressel & Condon (119781 , compiled in the NED database. The datum in the soft 
X-ray range come from the RO SAT All Sky Survey, reported by Walter & Fink ( 1993 1 and the 
observed flux at 1450 A comes from IUE observations of Rodrfguez-Pascual et al. J1997i . 

Although the observed SED was constructed from data obtained at very different epochs, it 
should be representative thanks to the limited variability displayed by Ark 564 (Romano et al. 
120021 . 
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Fig. 1. The observed continuum of Ark 564 compared to the results of the multi-cloud model 
Ml. The observed data from the NED are shown as filled squares. Open squares refer to the 
data of Walter & Fink (1993) in the X-ray and of Rodriguez-Pascual et al. dl997i in the UV. 
The multi-cloud emission, corresponding to dust emission and bremsstrahlung (thick solid lines) 
obtained from the weighted average of the clouds (thin solid lines) are shown, as well as the 
synchrotron emission from the shock front (dashed line). The old stellar population is represented 
by a blackbody (dot-dashed line). The frequency range where intrinsic absorption reduces the 
continuum flux is limited by the two thin vertical lines (see text). 

The most important physical mechanisms that should contribute to the emission of the con- 
tinuum in the different frequency ranges are (1) bremsstrahlung from the photoionized gas (v ~ 
5xl0 13 - 5xl0 15 Hz), (2) bremsstrahlung from the hot gas in the downstream region of shocked 
clouds (v > 10 16 Hz), (3) emission from the old star population (v ~ 5xl0 13 -10 15 Hz), (4) dust 
re-radiation in the infrared range (v ~ 10 12 -10 14 Hz), and (5) radio emission (v < 10 11 Hz). 

Besides the line and continuum spectra, the SUMA code calculates re-radiation by dust in 
a consistent way. Dust is heated by the nuclear radiation flux; moreover, gas and dust mutually 
heat each other by collisions in shock turbulent regimes (Viegas & Contini 1994 1. Therefore, the 
frequency corresponding to the dust re-radiation peak depends on the shock velocity. At high V s 
the grains will reach high temperatures in the post shock region and dust re-radiation will peak in 
the mid-IR, while at low velocities the peak appears in the far-IR. The peak intensity depends on 
the dust-to-gas ratio. Sputtering of the grains, which depends on the shock velocity, is calculated 
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throughout the clouds. Small grains are rapidly sputtered. Therefore, we adopt silicate grains 
with an initial radius of 0.2 /mi. 

Line and continuum spectra must be modeled consistently. The single-cloud models are cal- 
culated assuming a standard dust-to-grain, by number, d/g - 1CT 15 , providing a first fit to the 
emission lines. Once the best fit is reached with a multi-cloud model, the continuum emission is 
checked and the d/g values for the single-clouds may change. In this case, an iterative process is 
carried out, in order to obtain a self-consistent model accounting for the observed emission-line 
and continuum data. Notice that an increase of d/g in the postshock region, before complete sput- 
tering of the grains, acts as an increased density, speeding up the cooling processes. However, 
an increase by a factor less than 50 changes the cloud continuum emission without affecting the 
emission-line intensities. 

The relative weights adopted to calculate the continuum SED are the same as those used to 
model the line spectra. The dust-to-gas ratios (by number) that fit the dust re-radiation fluxes in 
the IR appear in Table 13 At high V, (models ml 3 and ml 4) the grains are completely sputtered 
at the shock front and we cannot determine the d/g ratio. The dust-to-gas ratio, by number, is 
<10~ 13 in low velocity clouds and as high as 8xl0~ 13 in clouds with V s =500 kms -1 . Recall 
that 10~ 14 by number corresponds to 4xl0~ 4 by mass adopting silicate grains with a density of 
3 gr cirT 3 . It means that in Ark 564 some high velocity clouds have dust-to-gas ratios higher 
than values found in the Galaxy. The high d/g found for some models (see Table [5J agrees 
to the relatively high E(B-V) considered by Crenshaw et al. (2002i. They claim, in fact, that 
the continuum spectrum of Ark 564 is clearly much flatter than that of Mrk 493 in the UV. This 
suggests that the spectrum of Ark 564 may experience a significant amount of internal reddening, 
which is supported by the Ha/H/3=3A measured by E97 and the small ratio of UV to X-ray 
continuum flux observed in Ark 564 compared to most Seyfert Is (normal and NLSls). 

The two thin vertical lines in Fig. ^ show the range of frequencies corresponding to strong 
absorption in the UV (see Contini et al. 2002a I, Within this interval, intrinsic absorption reduces 
the fluxes. Notice, for instance, that the hydrogen column density is > 10 22 crrT 2 for model ml4. 
The data at v < 10 15 Hz show that contribution from the old stellar population is present in 
Ark 564 (dash-dot line), as is found for most AGN. In the radio band, only one measurement in 
one frequency is available, so we cannot decide whether there is also some contribution from 
synchrotron emission created by Fermi mechanism at the shock front (dashed line, see Contini 
& Viegas DHS . 

5. Discussion 

The models which could better fit the coronal lines of NLS1 galaxies, including Ark 564, were 
already discussed in RVPP02 and they are based on line ratios extracted from a general grid of 
composite models (Contini & Viegas 2001 1. The asymmetries toward the blue observed in the 
line profiles of Ark 564 (E97) suggest that these lines originate in a gas outflowing from the 
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nucleus, supporting the models that combine the effect of shocks and photoionization. That first 
approach was a start for a more detailed analysis of the NLR of these objects, as presented in this 
paper. Here, emission-lines in a large range of wavelengths and ionization levels constrain the 
modeling and reveal that many different conditions coexist in the NLR. Notice that multi-cloud 
models were used previously in the detailed modeling of NGC5252, Circinus, NGC4151, and 
NGC7130 (Contini et al. 1 1 998al fl998bl I2()()2al I2002bl 

The multi-cloud models in Ark 564 show three main types of cloud velocities: (a) Those 
between 100 and 150 kms -1 , which are about the lowest observed in AGN and are typical of 
LINERS. Very low velocities clouds, with V s ~ 60 km s _1 , have been found to contribute only to 
the spectra of Circinus; (b) Velocities of about 300-500 kms -1 , generally present in the NLR; 
(c) High velocity clouds (> 1000 km s~'), also invoked in modeling the Seyfert galaxies referred 
above, have been used in order to explain high ionization lines ( e.g. [Fex]) as well as the soft 
X-ray emission. In the high V s clouds the d/g ratios are higher than in lower velocity clouds by a 
factor of ~ 100 both in Circinus and in NGC4151. High velocity shocks lead to a strong heating 
of gas and dust. However, the upper limit showed by Ark 564 in the NIR continuum at log v = 
13.5 constrains the d/g in the high velocity clouds to lower values. The same conclusion was 
obtained for NGC 5252. 

Preshock densities «o ^ 1000 cirT 3 are found in Ark 564, as well as in the other previous 
modeled galaxies. Only in the warm absorber of NGC 4051, another NLS1, preshock densities 
up to 10 6 cirT 3 were found in clouds characterized by very small D (> 2. 10 12 cm). 

Although the models for Ark 564 were chosen among those that better fit the strongest lines, 
they should be able to reproduce as much lines as possible. We briefly show the complex mod- 
eling of such a rich spectrum, discussing for instance [Si x]/[Si vi] relative to [Fe xi]/[Fe x] line 
ratios. The observations show that [Six]/[Sivi] > 2 and [Fexi]/[Fex] > 1. Model m6, which is 
matter-bounded and RD, gives the right trend for the Fe lines but strongly underestimates the 
[Si vi]/Hy6 line ratio, which is otherwise high in models composed of clouds with large D. In 
order to get a consistent picture, the models must be summed up with a fine tuning of the relative 
weights in order to reach an agreement between observed and calculated ratios for all the lines. 

Finally, if the NLR spectrum of Ark 564 is well reproduced by composite models accounting 
the effect of photoionization by a central source and shocks, it is then natural to ask about the 
origin of such shocks. Although SUMA does not discriminate about the shock sources, the fact 
that Ark 564 has an extended radio emission makes that jet-induced shocks be a likely scenario. 
This hypothesis is particularly suitable to explain the low-ionization lines. High velocity shocks, 
found to be adequate for coronal lines, may have their origin from outflowing material that evap- 
orates from the torus. This last hypothesis, in fact, was already considered by E97 and RVPP02 
from their observations of coronal lines. In summary, the broad range of shock velocities found 
from our modeling is consistent with shocks generated, at least, from two distinct mechanisms. 
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6. Concluding remarks 
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In the previous sections we have modeled, in detail, the observed continuum and line emission 
spectrum of the NLS1 galaxy Ark 564. Particular emphasis is put on the constrains imposed by 
the coronal lines. This is the first time that such an approach is carried out for the NLR of a NLS 1 
galaxy. The line spectrum ranges from UV to NIR wavelengths while the continuum extends 
from radio to X-ray. Composite models accounting for the combined effect of photoionization 
by a central source plus shocks are used to that purpose. We found a multi-cloud model that 
can suitably reproduce most of the observed line ratios within a factor of 2 and yields shock 
velocities consistent with the observed line widths. The picture that emerges from our modeling is 
a NLR composed of cloud showing a large range of velocities (V s = 100 - 1500 km s _1 ), preshock 
densities («o = 100 - 700 cirT 3 ) and illuminated by different flux intensities, from in SD models 
to Fh- 5x10 12 photons cirT 2 s _1 eV~' at 1 Ryd in RD models. The cloud geometrical thickness 
D was found to be between 3 and lOpc, but some of them represent filaments with D < 1 pc. The 
single-cloud models are summed up adopting relative weights. The best fitting model (Ml) is 
obtained by adopting the highest weights for SD clouds with V. s =150 and high velocity-density 
RD clouds reached by a strong radiation flux. These high velocity clouds, particularly SD clouds 
with V s = 1500 kms~' account for the soft X-ray emission and for most of the flux from the 
[Fe x], [Fe xi], and [Fe xm] lines (70%, 60%, and 55%, respectively). Clouds with V,=150 km s _1 
and «o=600 cirT 3 , illuminated by the strong central continuum flux, also contribute to the coronal 
lines. On the other hand, low ionization lines in the optical range are mostly emitted from SD 
clouds. In the modeling process we have considered forbidden line intensities relative to [Om], 
in order to avoid the broad line contribution that is present in the permitted lines. 

The input continuum energy distribution in the UV-X-ray region that photoionizes the NLR 
of Ark 564 is characterized by a spectral index a=-1.5. This value is very close to the a measured 
from UV and X-rays observations. That means that the NLR of Ark 564 is directly illuminated by 
the radiation of the central source. In addition, shocks, probably originated from the interaction 
between the radio-jet and NLR gas, provides an additional input mechanism for the production 
of coronal lines. 

To complete the modeling, we have roughly evaluated the contribution of the BLR to the total 
H/? line and to the permitted lines in the UV. A quantitative calculation is not possible because 
the SUMA code is not adapted to the calculation of lines from high density gas. However, it is 
found that Hfibmad/HPnarrow should range between 1 and 2, in very good agreement with the value 
found by Rodriguez- Ardila et al ( 2000 1 for NLS 1 galaxies and determined from Gaussian fitting 
to the observed permitted profiles. 

Finally, we could consistently explain the data of the continuum by the multi-cloud model 
Ml. A gas-to-dust ratio between 10~ 15 and 10~ 12 by number has been found for the NLR clouds, 
showing different values in different clouds. 
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Table 2. Single-cloud models 





ml 


m2 


m3 


m4 


m5 


m6 
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m8 


m9 


mlO 


mil 


ml2 


ml3 
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100 


100 
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150 


300 


500 
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50(1 
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30 
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300 


F H ' 
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1.(11) 


7.(11) 
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1.(12) 


3.(12) 


- 


1.(10) 


5.(12) 


- 


1.(11) 


5.(12) 


- 


D 4 


1. 


3. 


3. 


3. 


0.003 


0.01 


0.003 


3. 


3. 


3. 


3. 


1. 


0.1 


3.85(-2) 




1.4(-2) 


10(1) 


4.6(-l) 


6.1(-1)1 


1.6(-5) 


3.3(-3) 


3.7(-l)7 


1.8(-4) 


1.4(-3) 


2.49(+l) 


1.0(-2) 


8.7(+3) 


1.3(-2) 


4.5(-2) 


Model Ml : 






























wr 
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- 


- 


38. 


21. 


- 
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- 
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- 
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- 
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- 
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- 
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- 


- 


- 


8.6 


- 
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- 


- 
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- 
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- 


- 


- 
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- 
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- 
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- 
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%([Ci]) 


- 


- 


- 


- 
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- 


- 


- 
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- 
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- 


- 
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4. 


15. 
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11. 
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23. 






1.2 
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%([Sivi]) 
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11. 
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; 
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; 
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- 


- 
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- 


- 


- 


i. 


- 
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- 


- 
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19. 


- 
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- 


- 


- 
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- 
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1.1 
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- 
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- 
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- 


- 
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- 


- 
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- 
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- 
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- 


- 


16. 
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- 


- 
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- 


- 


- 
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22. 
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- 


- 
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- 


- 
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Model M2 : 






























wr 
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2.4(-4) 


1.6(-5) 


1.6(-4) 
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4.0(-2) 


1.2(-3) 


1.6(-1) 
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- 


- 


19. 


- 
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- 


- 
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- 
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- 
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- 
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- 


- 


- 
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- 
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- 


- 
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0.5 
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- 


- 
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- 
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9.3 


3.7 
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4.2 
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4. 
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1.5 
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3. 


5. 
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1.1 


0.7 


11. 






0.8 
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75. 


2.6 


3. 






%([Nev]) 








7. 




6.1 


69. 
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13. 
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Col. 1: lines 1-4 show the input parameters, followed by the H/3 calculated absolute flux for each model and then the output for the multi-cloud model Ml (relative 
weights wr, dust-to-gas ratios and percent contribution of each model to each line); the same for the multi-cloud model M2 appears in the last lines of the table. 
The djg results come from the SED which is modeled only for the selected model Ml. 
Cols. 2-15 : the output of single model predictions. 

1 in kms" 1 

2 ■ -1 
'incm J 

3 in photons cm" 2 s" ] eV" 1 at 1 Ryd 

4 in 10 19 cm 

5 in 1(T 15 

6 in erg cm -2 s -1 
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Table 3. Multi-cloud models results 
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[S in] ,19069+ 


14.0 


10.0 


19.0 


0.71 


1.36 


0.94 


0.86 


[Ci] 0.985//m 


0.19 


0.20 


0.45 


1.01 


2.5 


1.33 


1.5 


[Svm]0.991/im 


1.5 


1.7 


2.0 


1.18 


1.4 


1.55 


0.87 


[Fexm] 1.074/im 


2.4 


2.3 


5.0 


0.96 


2.0 


1.3 


1.3 


[S IX] 1.252//m 


1.7 


0.9 


0.7 


0.52 


0.43 


0.68 


0.27 


[Six] 1.430//m 


4.0 


2.0 


3.7 


0.56 


1.0 


0.74 


0.62 


[Si VI] 1.963/im 


1.6 


1.2 


1.9 


0.76 


1.2 


1.0 


0.76 


[Olll]////3 


1.19 


1.56 


0.75 










Ha/HP 


2.99 


2.77 


2.76 











1 RM1 and RIVE stand for the ratio {lxllOm\) calc l(l x l\Om\) obs for Ml and M2, respectively. 

2 RH/31 and RH/J2 stand for the ratio (lA/Hfi)calc/0-A/Hfi)<,b s iot Ml and M2, respectively. 
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Table 4. Comparison {hlWlcaicl{hlW) l obs for the UV lines. 



line Ml 



Lya 1.0 

Nvl214 0.33 
Oiv]+Siiv ,11400 1.27 

N iv] ,11486 0.50 

Civ ,11550 1.23 

Hen ,11640 4.7 

Om] ,11663 0.66 

Nm] ,11750 0.22 

Sim] ,11890 1.60 

Cm] ,11909 1.19 

[Ne iv] ,12423 1.36 

Mg n,12800 0.20 



1 from Crenshaw et al. ( 2002 Table 2, reddening corrected) 



